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ABSTRACT 

We report on the multi-wavelength observations of PKS 1510-089 (a flat spectrum 
radio quasar at z=0.361) during its high activity period between 2008 September and 
2009 June. During this 11 months period, the source was characterized by a complex 
variability at optical, UV and 7— ray bands, on time scales down to 6-12 hours. The 
brightest 7— ray isotropic luminosity, recorded on 2009 March 26, was ~ 2 x 10 48 erg s _1 . 
The spectrum in the Fermi-LAT energy range shows a mild curvature well described by a 
log-parabolic law, and can be understood as due to the Klein-Nishina effect. The 7— ray 
flux has a complex correlation with the other wavelengths. There is no correlation at all 
with the X-ray band, a weak correlation with the UV, and a significant correlation with 
the optical flux. The 7— ray flux seems to lead the optical one by about 13 days. From 
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the UV photometry we estimated a black hole mass of ~ 5.4 x 1O 8 M0, and an accretion 
rate of ~ 0.5 Mq/yt. Although the power in the thermal and non-thermal outputs 
is smaller compared to the very luminous and distant flat spectrum radio quasars, 
PKS 1510-089 exhibits a quite large Compton dominance and a prominent big blue 
bump (BBB) as observed in the most powerful 7— ray quasars. The BBB was still 
prominent during the historical maximum optical state in 2009 May, but the optical/UV 
spectral index was softer than in the quiescent state. This seems to indicate that the 
BBB was not completely dominated by the synchrotron emission during the highest 
optical state. We model the broadband spectrum assuming a leptonic scenario in which 
the inverse Compton emission is dominated by the scattering of soft photons produced 
externally to the jet. The resulting model-dependent jet energetic content is compatible 
with a scenario in which the jet is powered by the accretion disk, with a total efficiency 
within the Kerr black hole limit. 
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Introduction 



Among blazars, flat spectrum radio quasars (FSRQ) are those objects characterized by promi- 
nent emission lines in the optical spectra. The typical spectral energy distribution (SED) of blazars 
has a two bump shape. According to current models the low energy bump is interpreted as syn- 
chrotron emission from highly relativistic electrons, and the high energy bump is interpreted as 
inverse Compton (IC) emission. In FSRQs the IC bump can dominate over the synchrotron one 
by more than an order of magnitude. It is widely believed that in these sources th e IC compo 



nent is dominated by the scatteri ng of soft photons produced externally to the jet (jSikora et al 



1994 



Dermer &: Sch lickeise rl 120021) . rather th an by the synchrotron self Compton emission (SSC, 



Jones et al.lll974l : lGhisellini fc Maraschilll989i ), In the external radiation Compton scenario (ERC), 
the seed photons for the IC process are typically UV photons generated by the accretion disk 
surrounding the black hole, and reflected toward the jet by the broad line region (BLR) within a 
typical distance from the disk in the sub pc scale. If the emissio n occurs at larger distances, the 
external radiation is likely to be provided by a dusty torus (DT) (JSikora et al.ll2002l ). In this case 
the radiation is typically peaked at IR frequencies. 

The study of the SEDs of blazars and their complex variability has been greatly enriched since 



the 2 008 August start of scientific observations by the Large A rea Telescope (LAT, lAtwood et al 



2009) on the Fermi Gamma-ray Space Telescope (JRitzl 120071 ) . thanks to its high sensitivity and 
survey mode. 

One of the most active blazars observed in this period was the FSRQ PKS 1510-089. This 
object has an optical spectrum characterized by prominent emission lines overlying a blue continuum 



(JTadhunter et al.lll993l ) at a redshift z = 0.361 ([Thompson et al,l ll990). Radio images show a bright 
core with a jet which has a large misalignment bet ween the arcsecond and milliarcsecond scales; 
superluminal velocity up to ~ 20c are also reported (JHoman et al.l 12002 ). 



PKS 1510-089 was already detected in 7— rays by EGRET (JHartman et al.lll999i ) and exhibited 
a very interesting activity at all wavelengths. It was also detected by AG ILE during ten days of 
pointed observations from 2007-08-23 to 2007-09-01 (jPucella et al.l bo08h . In the period 2008- 
2009 PKS 1510-089 was observed to be bright and highly yariable in several freque ncy bands. 
In gamma rays it was detected in 2008 March by AGILE (JD'Ammando et al.l |2008| ) and other 
bright phases were observed in the subsequent months by both Fermi-LAT and AGILE (ITramacere 



2008 



Cutini k, Havs 



Villata et al. 



Ciprini fc Corbelll2009l : iD'Ammando et al.ll2009bl : IPucella et al.ll2009l : IVercellone et al 



2009). High states in X-rays and i n optical were reported by_ 



(2009a|), and lLarionov et al. 



a. 



In a recent paper, iMarscher et al. 
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presented data from a multi-wavelength (MW) campaign concerning the same flaring period of 
PKS 1510-089. In that paper, the authors focus on analysis of the parsec scale behavior and 
correlation of rotation of the optical polarization angle with the dramatic 7— ray activity. 

In the present paper we describe the results of the LAT monitoring together with the related 
multi-wavelength campaigns covering the entire electromagnetic spectrum. We present a detailed 
analysis of the 7-ray spectral shape and spectral evolution, and of the MW SED modelling and 
interpretation. This paper is organized as follows: in Section 2, we report results on the 7-ray 
observation of PKS 1510-089, and we study the spectral shape and its evolution. In Section 3, we 
summarize multifrequency data obtained through simultaneous optical-UV-X-ray Swift observa- 
tions, and radio-optical observatories. In Section 4 we present the results of the multifrequency 
data and their connection with the 7— ray activity. In Section 5 we report our conclusions about the 
MW data, and we use a phenomenological analysis to estimate some of the physical fundamental 
parameters, such as the black hole mass, the accretion disk bolometric luminosity, the shape of the 
electron distribution and the beaming factor. We then model the observed SEDs, and we comment 
on the jet energetics. Furthermore, we compare PKS 1510-089 with other powerful FSRQs observed 
by Fermi. In Section 6 our final remarks are reported. 

In the follo wing we use a ACDM (concordance) cosmology with values given within 1 a of the 
WMAP results (JKomatsu et allboogh . namely h = 0.71, tt m = 0.27, and Q A = 0.73, and a Hubble 
constant value Ho=WOh km s _1 Mpc -1 , the corresponding luminosity distance (djj is ~ 1.91 Gpc 
(~ 5.9 x 10 27 cm). 



2. Fermi-LAT data and results 



The LAT data presented here were collected from 2008-08-04 to 2009-07-01. Only events with 
energies greater than 200 MeV were selected to minimize the systematics uncertainties. To have 
the highest probability that collected events are photons the diffuse class selection was applied. A 



further selection on the zenith angle > 105° was applied to avoid contamination from limb gamma 
rays. The analysis was performed using the Science Tools package! (v9rl5p5). The Instrument 
Response Functions (IRF) P6_V3_DIFFUSE were used. These IRFs provide a correction for the 
pile-up effect. To produce light curves and spectral analysis the standard tool gtlike was used. The 
photons were extracted from a region of interest (ROI) centered on the source, within a radius of 7°. 
The gtlike model includes the PKS 1510-089 point source component, and all the point sources 



form the first LAT catalog (jAbdo et al.ll2010al ) that fall within 12° from the source. The model 
includes also a background component of the Galactic diffuse emission and an isotropic component, 
both of which are the standard models available from the Fermi Science Support CenteqH (FSSC). 
The isotropic component includes both the contribution from the extragalactic diffuse emission and 
from the residual charged particle backgrounds. The estimated systematic uncertainty of the flux 
is 10% at 100 MeV 5% at 500 MeV and 20% at 10 GeV. 



2.1. Temporal behavior 

We extracted light curves from the entire data set, to investigate the flaring activity. To take 
into account possible biases or systematics when the source flux is faint, we used two different time 
binnings of one day and one week. The light curves were extracted using gtlike , fitting the source 
spectrum by means of a power-law distribution (dN/dE oc E~ a ~>), where a 7 is the photon index, 
following the prescription given in the previous Section. The flux was evaluated by integrating the 
fitted model above 0.2 GeV. The lower panel of Fig. Q] shows clearly four major flaring episodes: 
between 2008-08-30 and 2008-09-26 (flare a), between 2009-01-04 and 2009-01-27 (flare 6), between 
2009-03-10 and 2009-04-09 (flare c), and between 2009-04-15 and 2009-05-12 (flare d). The source 
was almost quiescent between the end of 2008 Sep. and the beginning of 2009 Jan. (see. Tab. [I] for 
a summary). The flux light curves with different temporal binning are compatible, with the daily 
integration binning showing better rapid flux variations that are smoothed in the weekly binning. 
A study of t hese variations based on the autocorrelation, Fourier analysis and structure function is 



presented by lAbdo et al.1 (j2010bl ) together with other blazars. 



Fig. [2] shows close-up light curves of the four flares and the green line represents the optical 
data in the R filter (see Sect. [3]and[I|). Since the statistics during the flares were high, it was 
possible to use also a 12 hr binning (blue points). Typically, the flares have a complex structure 
with peaks having durations from about 1 to 5 days and their moderately asymmetric profile can 
result from the overlapping of subsequent episodes. In a few cases, significant variations by a factor 
of 2 within 12 hr were detected. To have a better estimate of the rising and decaying timescales, 
we fit two rapid flares, with an almost regular shape, using an analytical law of the form A ■ exp*' T . 
In the case of the flare peaking at t ~ 54846 MJD (dashed black line in panel b of figure [2]) we find 



* http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
' http://fermi.gsfc. nasa.gov/ssc/data/access/lat/BackgroundModels. html 



a rise faster than the decay: the rising timescale is r ~ 0.3 days and the decaying one is ~ 1.4 days. 
The second event (flare peaking at t — 54962 MJD, dashed black line in the d panel of [2]) followed 
the opposite behavior, having a best fit rising timescale of ~ 0.8 days, and a decay time of ~ 0.25 
days. Using a bin width of 6 hr, the shape of the flare is nearly symmetric, with rise and decay 
e-folding time of about ~ 0.12 days. Such a fast variability can constrain the radiative region size 
Rrad by the well known relation 



J^rad — 



c At S 

TT7' 



(i) 



where c is the speed of the light, 5 = 1/(T(1 — (3 cos 6)) is the beaming factor depending on the bulk 
Lorentz factor T and a viewing angle 9 ~ 1/I\ and z is the c osmological redshift. Adopting the 
very fast superluminal velocity reported bv lHoman et al.l (J2002J) ((3 app ~ 20, in very good agreement 
with the results presented in Section I4.3P from 5 ~ T ~ j3 app we can estimate R ra d i$ 9 x 10 15 cm 
(in the case of r ~ 0.25 days). We will compare this result with other constraints derived in Sect. 
EH 



2.2. Gamma-ray Spectra 



We analysed the 7-ray spectral shape of PKS 1510-089 during the whole period, the quiescent 
state, and the four flaring episodes using three spectral models: a power-law (PL) a log-parabola 
(LP), dN/dE oc #/£- a -r-/31og(£/Eo) (JLandau et al.lll986l : lMassaro et alJl2004J ). and a broken power- 
law (BPL). In the case of LP spectral law, the parameter f3 measures the curvature around the 
peak. The LP distribution has only three free parameters, and the choice of the reference energy 
Eq does not affect the spectral shape; we fixed its value to 300 MeV. We performed the spectral 
analysis using an unbinned maximum-likelihood esti mator (gtlike) and the same prescription 
given in Sec. We used a likelihood ratio test El (LRT. lMattox et al.lll99a ) to check the PL model 
(null hypothesis) against the LP model (alternative hypothesis). Since the PL is often rejected, 
we test also the LP model (null hypothesis) against the BPL model (alternative hypothesis). The 
results concerning the LRT are summarized in Tab. [2j and in Tab. [3] we report the details of the 
spectral analysis for each time range and spectral model. Due to the non derivable character of 
the BPL law, we used also the loglikelihood profile method to determine the best fit parameter 
for this model. The corresponding statistical uncertainty was estimated from the difference in the 
likelihood value w.r.t. its minimum such that —2AL=1. The 7-ray spectrum of PKS 1510-089 is 
well described by a log-parabola, with the only exception being the quiescent state. The value of 
the LRT reported in Tab. [2] shows that both for the flares (a, b, c, d) and for the full period, the 
LP model describes the spectrum better than PL with a probability higher than ~ 99.6%. The 



'The LRT statistic is defined as: LRT = -2Log(L0/Ll), where L0 and LI are the maximum likelihood estimated 
for the null and alternative hypothesis respectively. 



BPL, on the contrary, does not provide an improvement w.r.t. the LP model. The only exception 
is the flare b, but the false positive probability is about 27.33%, so there is no evidence to reject 
the null hypothesis. The LP model is then preferred because of the lower number of parameters. 
Moreover the curvature parameter (3 can be linked to physical processes such as the acceleration or 
the effects of the Klein-Nishina regime in the IC process, as we will discuss in Sect. I5.ll and l5.2i For 
a better visualization of the SED shape and to show the departure from a PL trend, we produced 
an SED by performing an independent likelihood analysis starting form a grid of 20 energy bins 
logarithmically equispaced. The bins were then grouped in order to have at least 10 photons per 
bin, and the highest energy bin was chosen according to the maximum energy encircled within 95% 
of the PSF. The results are shown in Fig. [3] and Fig [31 In Fig. [3] we show the full-period SED and 
the one extracted during the quiescent state. We plot by a dotted line the PL model, by a dashed 
line the LP model, and by a dot-dashed line the BPL model. With a red upward arrow we indicate 
the highest energy event within 95% of the PSF for the whole period data set, corresponding to 
an energy of approximately 30 GeV. From the plot of the PL model residuals (lower panel) it 
is possible to see clearly the departure from a PL trend: the deviations, both at low and high 
energies, suggests for a spectral curvature confirming the LRT results. Moreover, it is possible to 
note that the BPL model does not deviate significantly from the LP trend, supporting again the 
LRT analysis. The SEDs of individual flares are plotted in Fig. [Hand show that the spectrum was 
curved also during the single flaring episodes. We note that the flare-integrated spectral shape did 
not change significantly, despite the huge flux variations. 



2.2.1. Spectral evolution. 

To complete the analysis of the spectral behavior of PKS 1510-089, we investigate whether 
spectral changes are seen between the quiescent and the flaring state. Since we are mainly interested 
in the search of possible trends rather than in the best description of the spectral distribution, we 
simply evaluated the PL spectral indices in the various brightness states, which can be considered 
representative of the mean slope. In the upper panel of Fig. [5] we plot the photon index against the 
flux above 0.2 GeV, resulting from the same spectral analysis showed in Fig. [H In the case of daily 
integration (green circles), and more marginally for the weekly integration (red circles), this plot 
seems to show a softer when brighter trend, up to a flux level of F(E > 0.2GeV) ~ 2.4x10 ph 
cm~ 2 s _1 . Above this value of the flux the source has a harder when brighter trend. We analyse the 
correlation of the harder when brighter trend, for the weekly binning, using a Monte Carlo method 
that takes into account the dispersion of flux and index measurements. In detail, we re-sample 
the flux and index values for each observed pair, extracting the data from a normal distribution 
centered on the observed value, and with the standard deviation equal to the 1 a error estimate. 
We find a correlation coefficient of r = 0.43 with a 95% confidence limit 0.24 < r < 0.58. The trend 
for F(E > 0.2GeV) > 2.4xl0~ 7 ph c m" 2 s~ 1 is reported i n the inset of the upper panel of Fig. [5j 



Although for some EGRET blazars iNandikotkur et al.l (|2007l ) observed a similar flux-hardness 



anticorrelation at low fluxes, we need to take into account possible effects coming from the poor 
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statistics when the source flux is low. As first, we note that moving to the weekly binning the 
trend is less evident, although the flux range is the same as that of the daily binning. As a further 
check, in the lower panel of the Fig. [5] we plot the photon index against the number of photons 
predicted by the best-fit model. It is clear that the dispersion of the photon index is related closely 
to the number of predicted photons, above N ~ 20 the trend is the same for both the 2 integration 
timescales, and the photon index clusters around 2.5, without showing the very soft and very hard 
indices values present at low number of events. 

In conclusion, we can not exclude the presence of a softer when brighter trend, for low flux 
levels and short timescales, but the statistical effects we present do not allow to obt ain a purely 



physic al interpretation. Similar results have been found for other Fermi Blazars (jAbdo et al 



2010dl ). independently of their redshift or class (BL Lacs/FSRQs). 



Even if we did not find a strong evidence for a index-flux correlation, the variation of the 
photon index returned by the gtlike fit as a function of time (see upper panel of FigQ]), shows 
that the dispersion on the photon index is larger until mid March roughly, and gets narrower after. 
This feature is emphasized in Fig. [6J where in the upper panel we plot the histogram of the photon 
index for the weekly integration, before MJD 54905 (corresponding to 2009-03-15, blue shaded 
histogram), and after MJD 54905 (red empty histogram). In the lower panel we plot the same 
analysis for the case of daily integration. The distributions before and after MJD 54905 have the 
same mean, in both daily and weekly integrations (~ 2.5), but very different standard deviations. 
In the case of daily integration we have 0.45 and 0.21 before and after MJD 54905, respectively. 
In the case of the weekly integration we have a standard deviation of 0.28 before MJD 54905 and 
0.07 after. To test more quantitatively whether or not the distributions are different, we use a 
Kolmogorov-Smirnov (KS) test. We applied the test to the distributions of the spectral indices 
before and after MJD 54905, for the daily and weekly integrations. The test returns a p- value of 
~ 0.13 and ~ 0.35, for the case of weekly and daily binning respectively. The KS test gives only a 
marginal indication that the data sets, before and after MJD 54905, may not be drawn from the 
same distribution. 

In conclusion, the typical 7-ray photon index of PKS 1510-089 is around 2.5 and values largely 
different from this were never observed in high states. LP best fits indicate a significant but quite 
mild spectral curvature, that in the brightest flares (b, c and d) was always very close to /3 = 0.1. 



3. Multifrequency observations and data reduction 

The unique, high quality data provided by the LAT instrument can not be physically fully 
understood without simultaneous multifrequency observations. The spectral curvature and the 
spectral evolution observed in the 7-ray band, need to be compared to SED evolution from the 
radio to the hard X-ray. Radio data, and in particular VLBI data allow us to constrain the beaming 
factor and to cross check this result with that obtained from the 7— ray transparency. X-ray data 
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can shed light on the balance between the SSC and the ERC component, and allow us to estimate 
the spectral shape of the low energy branch of the electron energy distribution. UV data provide 
information about the Big Blue Bump (BBB) radiation, and combined with optical data tells us 
about the high energy branch of the electron distribution. Moreover, UV/optical data constrain 
the peak flux and energy of the low energy bump, determining the ratio between the output of the 
synchrotron distribution to that of the IC one. 

In the following of this subsections we report the reduction of the data collected at different 
wavelengths as a result of pre-planned campaigns (GASP optical-to-radio and VLBI radio data) or 
as ToO triggered by the LAT flaring activity (Swift data). In the next Section (Sec. HJ we discuss 
the multi-wavelength results and their connection to the LAT data. 



3.1. SWIFT-SAT and XRT data 



We analyzed XRT (IBurrows et al]l2005l : lGehrels et al.ll2004l ) data using the xrtpipeline tool 



provided by the HEADAS v6 . 7 software package, for data observed in photon counting mode. Events 
in the 0.3-10 keV energy band were extracted, selecting grades in the range 0-12, and default 
screening parameters to produce level 2 cleaned event files were applied. Due to the low count rate 
(< 2 cts/s) we did not find any signature of pile- up effect. 

We used data from the Burst Alert Telescope (BAT) on board the Swift mission to derive 
the spectrum of PKS 1510-089 in the 14-195 keV band. The spectrum is constructed by averaging 
the spectra of the source extracted over short exposures (e.g. 300 s) and it is representative of 
the sources emission over the 5-year time range spanned b y the observ a tions. These spectra are 
accurate to the mCrab level and the reader is referred to lAiello et al.l ( 20081 . l2009al lbl) for 



more 



details. PKS 1510-089 is bright in BAT and the approximate significance of the BAT spectrum 
used for this analysis is ~13<r. 



3.2. Swift-UVOT data 

We followed the steps outlined in the UVOT User's Guide, to perform UVOT data reduction 
and analyses in all the 6 available filters (V, B, U, UVW1, UVM2, UVW2). We started from the raw 
data stored in the HEASARC archive and we made sure that the sky coordinates were updated, the 
modulo-8 correction was applied, duplicated FITS extensions have been removed, and the aspect 
correction was calculated. Based on the AGN intensity, the optimal source extraction region is a 5" 
circle. The background region is an annulus with inner-outer radii of 15"— 27", 27"-35", depending 
on the filter used. In order to improve the astrometry, the NED position has been adjusted using 
the uvotcentroid task, and field of view sources have been excluded from the background region. 
The standard output of the uvotsource task has been used to extract the photometric light curves. 
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We corrected the magnitu des for Galacti c extin ction assuming E(B — V)q s ,\ = 0.097 mag. This 
value was calculated from ISchlegel et al.l (|1998l ) tables using tools provided by the NASA/IPAC 
archivq§. The absor ption for the other filters was calculated according to the extinction laws of 
Cardelli et al.l (J1989). The de-redde ned magnitu d es we re converted into fluxes in physical units 
taking into account the zeropoints by iPoole et al.l (|2008l ). 



3.3. Optical Near-IR and radio observations by the GASP 

The GLAST- AGILE Support Program (GASP) is performing a long-term monitoring of 28 7- 
ray loud blazars in the optical, near- infrared, mm, and radio bands (jVillata et al.ll2008l . l2009bl ). The 
GASP has been following PKS 1510-089 sin ce 2007 January, and contributed to multi- wav elength 
studies involving 7-ray data from AGILE (jPucella et al.l 120081 ; iD'Ammando et al.l l2009al ). The 
optical and near-infrared GASP data for the present paper were acquired at the following obser- 
vatories: Abastumani, Armenzano, Calar Alto, Campo Imperatore, Castelgrande, Crimean, Kitt 
Peak (MDM), L'Ampolla, Lowell (Perkins), Lulin, Roque de los Muchachos (KVA and Liverpool), 
Sabadell, San Pedro Martir, St. Petersburg, Talmassons, and Valle d'Aosta. Magnitude calibration 
was performed with respect to a common choice of reference stars in the field of the source from 
the photometric sequence by iRaiteri et al.1 (J1998I ) . Conversion of magnitudes into de -reddened flux 



densit ies was obtained by adopting the Galactic absorption value Ab = 0.4 16 from 



Schlei 



( 19981 ). consistent with the E(B — V) color excess, the extinction laws by ICardelli et al.l (1989), 
and the mag-flux calibrations bv lBessell et al.l (|1998). 



et al. 



The GASP mm-radio data were taken at Medicina (5, 8, and 22 GHz), Metsahovi (37 GHz), 
Noto (43 GHz), SMA (230 GHz), and UMRAO (4.8, 8.0, and 14.5 GHz). 



3.4. VLBI data 



The 2 cm VLBA / MO JAVE program ([Lister et al.l (|2009bl ) and references therein) has been 
monitoring PKS 1510-089 at 15 GHz with the Very Long Baselin e Array (VLBA) si nce 1995. 
Method of observations, data processing and imaging is discussed bv lLister et al.l (|2009bl ). Typical 
resolution of these images is about or better than 3 pc. 

In addition to the 15 GHz MOJAVE VLBA monitoring, single-epoch simultaneous multi- 
frequency 5-43 GHz VLBA measurements w ere done on 2009 April 09, in support of the first year 
Fermi observations (jSokolovsky et al.ll2010l ). Accuracy of flux density measurements is dominated 
by calibration uncertainties: about or less than 5% at 5, 8, and 15 GHz, about our less than 10% 
at 24 and 43 GHz. 



''The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under contract with the National Aeronautics and Space Administration. 
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4. Multifrequency Results and connection with the LAT data 

4.1. X-ray and Hard X-ray data 

We performed the spectral analysis of Swift-XRT data after grouping the photons to have a 
minimum number of 10 photons per bin, and we fitted the spectra by means of a photon power- 
law distribution F(E) = K E~ ax , plus a Galactic absorption with an equivalent column density 



Nh = 7.88 x 10 20 cm 2 ( Lockman fc Savage! 1 19951 ). We extracted the X-ray spectrum for each 



pointing. The corresponding spectral analysis results are reported in Tab. HI 

The X-ray light curve obtained from the fluxes reported in Tab. [5] shows a modest variability 
if compared to optical and 7-ray flares (see Fig. [7]). We prefer to present light curves in terms of 
vF(v) to make easier the comparison between the various bands and the SED changes. The average 
flux integrated in the 0.3 — 10.0 keV range is around 10 -11 erg cm -2 s" 1 , the lowest flux, recorded 
on 2009-01-16, was (6.5±0.8) xl0~ 12 erg cm" 2 s -1 . The highest flux, recorded on 2009-04-28, was 
(15.0±1.5)xl0 2 erg cm -2 s . In this case the flux increased by a factor of two within a day, 
and the spectrum reached the hardest state {ax = 1.13 ± 0.13). This is the most relevant X-ray 
flaring episode in our data set and looking at the MW light curve in Fig. [7] it seems to have no 
counterparts in other wavelengt hs. Since th e statistics are low, we performed the spectral analysis 



using the Cash method (C-stat) (|Cashlll979l ) based on the use of a likelihood function. This method 
returns flux and photon index values that are compatible with those coming from the x 2 method. 
Even if the two methods results are compatible, the significance of this flare is low (~ 2a), so we 
do not investigate possible physical implications. 

During our observations the source spectrum was always hard, with a photon index ranging 
between about 1.3 and 1.6. The plot of the flux in 0.3 — 10.0 keV range vs. the photon index 
(see Fig. |HJ) is compatible with a harder when brighter trend. Using the Monte Carlo method 
described in Sect. 12.2.11 we get a correlation coefficient r = —0.31 with a 95% confidence limit 
of —0.55 < r <__— 0.05. This spectral tr end is consistent with the same analysis performed by 



Kataoka et al.l ((2008). We note also that iKataoka et al.l (]2008l ) found a soft X-ray excess in the 



Suzaku data, but the statistics of the individual pointings in our data set are not sufficient to 
detect such a feature. 

In order to increase the statistics and to look for differences between the different 7— ray flares, 
we produced X-ray SEDs averaged during the b, c, and d 7-ray flaring intervals, and during the 
post-d flaring period, reported in Fig. [9) These SEDs show that the average state of the X- 
ray emission was almost steady, without drastic differences between the flares and the post-flare 
integration period. We also note that a possible soft X-ray excess is visible in the post- b flare 
averaged SED. The scatter plot in Fig. [10] shows no correlation between the XRT flux and the LAT 
flux. This absence of correlation is relevant to the understanding of the emission scenario that we 
will discuss in the Sect. [5j 

The 5-year integrated BAT SED is plotted in Fig. [9) The photon index, in the 14-150 keV 
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band, is 1.37^ 0ig . Despite the long integration time of the BAT data, the photon index value is 
almost compatible with the range of values observed in the XRT data in our data set, and in other 
historical observations. This suggests that the X-ray and hard X-ray flux and spectral shape of 
this source is quite stable, or at least that our X-ray sampling is representative of the X-ray and 
hard X-ray shape on timescales of years. 



4.2. Optical/near IR and UV data results 

Simultaneous Swift-UVOT and GASP observations provide a valuable data set to study the 
low energy bump of the SED. In Fig. QT] we plot the SEDs obtained from UVOT and GASP data, 
simultaneous within a daily timescale. These data show a minimum around the frequency of 
5 x 10 14 Hz, which does not seem to vary when the source is flaring. The frequency of the high 
energy peak can be well estimated and is close to 10 15 Hz, while that of the peak at low frequency 
cannot be well established and should be roughly estimated around 10 13 Hz. The UV peak, as in 
many quasars, is likely due to the the BBB that usually is understood as thermal emission from the 
accretion disk surrounding the black hole. Assuming that the disk luminosity is steady or slowly 
variable compared to the synchrotron emission, we expect that the UV excess gets less and less 
evident, as the optical flux increases. To test this scenario, in Fig. [12] we plot the ratio of the flux 
in the highest energy UVOT filter ([TVW2) to the flux in the B UVOT filter, as a function of the 
optical R flux. This plot shows that the UV spectrum gets harder when the optical R flux is lower. 
This trend is consistent with the UV excess decreasing as the R flux is increasing, as expected in 
a scenario in which the UV bump originates in the thermal emission of the accretion disk. Despite 
this trend, we note that even in the highest state plotted in Fig. [TT| the corresponding UV bump is 
still pr ominent. For other FSRQs, such as 3C 4 54.3, the BBB is not visible during the high optical 



states ( Giommi et al.1 12009 : IVillata et al.1 120061 ) . since it is dominated by the synchrotron flux, as 



we would expect here in the case of PKS 1510-089. In the panels a and b of Fig. [13] UV light 
curves for all the 6 filters are shown, and in the panel c we plot the ratio of the UVOT V filter 
flux to the UVOT W2 filter. The trend of the hardness ratio is clearly anticorrelated with that of 
the fluxes. Again, this supports the BBB scenario, with the UV spectrum harder when the BBB is 
more evident, namely when the synchrotron flux is lower. 

In Fig. Hwe report the scatter plot of the Fermi-LAT flux F(E> 200 MeV) vs. the UVOT 
i>F(v) in the UVVFi? filter, to check for a possible significant correlation between the 7— ray and UV 
fluxes. The correlation coefficient of the logarithms of the UV and 7-ray fluxes, obtained through 
the Monte Carlo method described in Sect. 12.2.11 is r = 0.2 with a 95% confidence interval of 
0.05 < r < 0.34. If we exclude the point with the lowest 7— ray flux, the correlation coefficient is 
r = 0.05, suggesting that the overall correlation is not significant. 

This lack of correlation hints that the spectral evolution of the UV spectrum results from the 
contamination of the high energy branch of the synchrotron emission, and that a change in the 
BBB luminosity is not the driver of the 7— ray luminosity variations. We will discuss this in Sect. 
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5. 

The optical (R band) and near-IR (J, H, and K bands) observations span the period from 
2008 May until the end of 2009 June. In Fig. [15] we report the light curves of the optical data. 
The best sampled band is the R one, with more than 600 observations. The optical flaring activity 
increased dramatically after mid March 2009, corresponding to th e 7-ray flare c. On 20 09-05-08 



the source was very bright with an R magnitude of 13.60 ± 0.02 (JLarionov et al.ll20 09a). and in 



two days it reached its historical peak at R = 1 3.07 zfc 0.02, about 1.3 magnitude brighter than 



the previous record level, observed on March 27 (JLarionov et al.ll2009bl ) of the same year. This 
dense monitoring is fundamental to understand the correlation between the optical and the 7— ray 
activity, both for the long-term trends and the single flares. In Fig. [16] we show the scatter plot 
of the Fermi-LAT flux F(E> 200 MeV) vs. the optical v¥(u) in the R filter. The correlation 
coefficient of the logarithm of the optical and 7-ray fluxes, evaluated through the Monte Carlo 
method described in Sect. 12.2.11 is r = 0.42 with a 95% confidence interval 0.36 < r < 0.46, higher 
than that found for the UV band. This finding hints that the driver of the flaring activity can be 
a change in the high energy branch of the electron distribution. We will investigate this scenario 
more accurately in Sect. [5] Despite the statistical significance of this result, there is an evident 
dispersion in the scatter plot, that could be related to the inter-band time lags. Indeed, temporal 
lags could be related to the internal source photon absorption, to the cooling time of the radiating 
particles, or to inhomogeneities in the emitting region. To search for possible time lag between 
the optical (R band) and the 7-ray band, we used the discrete cross correlation function (DCCF) 
method. Taking into account the whole data set, the DCCF analysis returns a lag of 13 ± 1 days, 
with the 7-ray band leading the optical one (see Fig. (XT] top panels). Since the R flux increased 
dramatically after the end of 2009 February, we divide the data set in two Sections, before and after 
MJD 54890 (corresponding to 2009-02-28). The DCCF results, plotted in the middle panels of Fig. 
\TT\ show that the 13 days delay is still apparent. The same result holds if we analyse individually 
the flaring sequences d and b (Fig. [TT\ bottom panels). Looking at the flare light curves it is clear 
that the 13 day lag is caused by a change in the relative flux of the two bands such that the 7-ray 
flux is stronger in the first part of the flare and the R-band flux gets stronger in the second half. 
We plot in the lower panel of Fig. [18] the LAT light curve above 200 MeV, and in the upper panel 
the optical light curve, backward shifted by 13 days: the four 7-ray flares a, b, c, d, seem to have 
an optical counterpart. In particular, the peaks of the bright events in May would be very close in 
time to those in the LAT daily light curve at MJD ~ '948 (see Fig. [2]). To check the effect of this 
lag, we evaluate the correlation between the logarithms of the optical and 7-ray flux after applying 
the above time shift (see Fig. [T9|) and find r = 0.62 with a 95% confidence interval 0.57 < r < 0.66. 
The values of the correlation coefficient and of its 95% confidence interval, after applying the 13 
days time shift, are significantly higher than those obtained without the time shift. Moreover, Fig. 
[19] and Fig. [TBI show that time shift reduces the dispersion in the scatter plot. 

This 13 day correlation between the optical and 7-ray emission appears in the different out- 
bursts during this flaring activity but we do not have any indication that it is a characteristic 
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behavior related to the temporal and energetic evolution of the flares. Only through a long-term 
simultaneous 7— ray and optical observations may we understand better the actual level of random- 
ness of this correlation and its possible physical meaning. 



4.3. Radio results 



Radio data have a much lower sampling w.r.t. optical and 7-ray, but at 14.5 GHz and 43 GHz 
it is possible to follow the overall trend of the radio flaring activity. Even if radio fluxes do not show 
a correlation with the optical and the Fermi-LAT flaring trend, it is possible to identify a large 
radio flare, starting at about MJD 54900 (see Fig. [7]). This radio flare, visible at 14.5, 43, and 230 
GHz, seems to starts quasi-simultaneously with the 7-ray flare c at 43 GHz, and keeps increasing 
until the end of flare d. The 230 GHz light curve shows a structure similar to a plateau, starting 
when the 7-ray flare c has ended, and a possible plateau is present also at 43 GHz, starting at the 
end of flare d. The 14.5 GHz fluxes seem to lag behind the 43 and 230 GHz. Analysis in the radio 
band is much more complex than optical and 7-ray data, because of the synchrotron self-absorption 
and the longer cooling times, and the possibility of the overlapping of different flares. 

In 2009 June-July the integrated parsec-scale flux density had reached its historical maximum 
since 1995 (~ 4 Jy), as the 15 GHz MOJAVE VLBA measurements show. The VLBI core flux 
density has also shown the highest value (see Fig. [20l upper panel). If this major radio flare, 
observed in the core of the parsec-scale jet, is connected to the huge 7-ray flares which have 
happened in the first half of 2009, this determines the source of the high energy emission to be 
located around the base of the parsec-scale jet following causality arguments, as suggested by 
Kovalev et al.l (120091 ). The delay between the peaks of light curves in the 7-ray and radio bands can 
be, at least partly, explained by the synchrotron self-absorption of the emission at radio frequencies. 

In Fig. [21] we present the three highest frequencies Stokes I parsec-scale images, from the 
5-43 GHz VLBA measurements performed on 2009 April 9. The size of the bright parsec-scale core 
at 24 and 43 GHz is estimated to be about 60-70 ^as or 0.3-0.4 pc. The core shows a flat radio 
spectrum (see Fig. 1221) indicative of a synchrotron self-absorbed region, while the first well resolved 
jet feature is already optically thin radio spectral index s r = —0.9 (F(u) oc u Sr ). 

The hi ghest apparent spee d of a component motion in the jet of PKS 1510-089 observed at 
15 GHz by iLister et al.l ()2009d ) is v app = 24 c which makes its j et highly relativist i c and Doppler 



boosted, a typical case for 7-bright blazars (JLister et al. 



2009 



let 

3). 



Marscher et al.l (|2010al ) report 



on two new knots, observed in VLBA images at 43 GHz. The first with an apparent speed of 
24 ± 2c, passed the core on MJD 54674.5 ± 20 (2008-07-27). The second knot passed the core on 
MJD 54958.5 ± 4 (2009-05-07) which is consistent, within the estimated uncertainties, with the 
huge optical flare observed on 2009-05-08. A feature, coincident to the second knot, is also seen 
to emerge in the MOJAVE 15 GHz VLBA images from 2009 June through 2009 Decembe r. The 
fitted speed at 15 GHz (675/uas y _1 ) is somewhat slower than the iMarscher et al.1 (|2010al ) speed 
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(970 ± 60 /ias y -1 ), but this could be due to blending with a third feature seen in the 43 GHz 
data to emerge sometime in June 2009 (JMarscher et al,ll2010bl ). Despite the high resolution of the 
43 GHz images, a more accurate ejection date for the latter feature could not be obtained due its 
rapid angular evolution and close proximity (0.1 mas) to the feature ejected in 2009 May. 



5. Interpretation and discussion 

5.1. Physical interpretation of 7-ray and MW data 

PKS 1510-089 was one of the brightest and most active blazars observed by Fermi-LAT during 
the first year of survey. Flaring episodes with timescales from weeks to months in addition to rapid 
and intense outbursts, were observed both at optical and 7— ray energies. The estimated isotropic 
luminosity above 100 MeV during the flare c was about 8xl0 47 erg s _1 . This value represents the 
brightest flare-averaged state of the source. The brightest daily time-resolved luminosity, recorded 
during the flare c on 2009-03-26, was of ~ 2 x 10 48 erg s . During the flare d, variability timescales, 
down to a fraction of a day, were observed both at optical and LAT energies. 



The rapid variability and the powerful 7— ray luminosity raise the problem of the pair pro- 
duction opacity. Indeed, without beaming effects, the source size estimated from the observed 
variability timescale (R ra d = cAt/(l+z)), make the source opaque to the photon-photon pair pro- 
duction process, provided that 7— ray and X-ray photons are produced cospatially. 

If the component emerged on 2009-05-07, observed also in the MOJAVE 15 GHz VLBA, is 
related to the flares a and b, then we can use the value of the beaming factor derived from the 
motion of radio knots, namely 5 ~ T ~ f3 app ~ 21. This estimate is compatible with other VLBI 
estimates (JHoman et al.1 |2002|) , and is slightly larg er than an alternat ive estimate based on the 
variability observed at 22 and 37 GHz, 5 var = 16.7 (JHovatta et al.ll2009l ). Taking into account the 
most rapid timescale estimated in Sect. 12. H At ~ 0.25 days, the actual emitting region size results 
of the order of R ra d < SAtc/(l + z) ~ 1 x 10 16 cm. Using the above fastest timescale and the quasi- 
simultaneous observed X-ray flux Fx — 8 x 10 _12 erg cm -2 s _1 (observed at a typical frequency 
of 10 18 Hz), one can impose a limit on the minimum value o f the beaming factor resulting in a 
source transparent to the p hoton-photon annihilation process (JMaraschi et al.lll992l : lMattox et al. 
1993J; iMadejski et al.l Il996l ) . Combining the source size from Eq. [TJ the X-ray photon energy in 
the source frame, and the intrinsic X-ray luminosity with the optical depth expression, we get a 
model independent estimate of 5 > 8. This lower estimate indicates that the values derived from 
the VLBI images are well above the pair production transparency limit. 

X-ray fluxes, spectral indices and trends in our data set are compatible with those reported i n 
previous analysis with several X-ray telescopes (BeppoSAX, ASCA, Suzaku: iKataoka et al.l l2008). 
The harder when brighter trend and the lack of X-ray /'y— ray correlation are useful to constrain 
both the low energy tail of the electron distribution and the emission scenario. The typical value 
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of the soft-to-hard-X-ray photon index ax is close t o 1.4, with a quit e narrow dispersion (~ 0.1 



and is similar to that observed in FSRQs with z > 2 (jPage et al.ll2005l ). Since in FSRQ objects the 
X-ray band samples the low-energy tail of the ERC component, the X-ray energy spectral index 
(sx) constrains t he slope p of the low-energy tail of the electron distribution in the range 1.6-2.0 



(j) = 2sx — 1, e.g. iRvbicki Sz Lightmanlll979l ) 



The high energy spectral index of the electron distribution can be estimated from the Opti- 
cal/UV spectral shape. In Sect. I4.2l we showed that the Optical/UV spectral index depends on the 
relative contamination between the synchrotron component and the BBB emission. A reliable esti- 
mate of the synchrotron spectral index could be achieved from the optical data, but, unfortunately, 
our spectral coverage with two or more simultaneous frequencies is only in a limited time window. 
In Fig. [23] we plot the spectral index in the R to H bandpass vs. the 7— ray flux integrated above 
200 MeV. If we take into account the blue points which refer to the subflare peaking at about MJD 
54909 (see figJ21 panel d), we note that the 7— ray flux increased by about a factor of 2.5 with 
the optical spectral index almost constant around ~ 1. This value of the optical spectral index 
corresponds to an electron energy distribution index of about 3. We note that the steadiness of the 
optical R — H spectral index, during this 7— ray subflare, is consistent with the spectral evolution 
of the 7— ray emission. Indeed, extracting the 7— ray photon indices simultaneous within one day 
with the optical ones, we note that also these values are almost stable around the value of 2.4. The 
decay of the flare occurred with a significant spectral softening, consistent with a cooling dominated 
regime, with a corresponding electron spectral index in the range between 3.5 and 4.0. 

The observed multi-wavelength SEDs of PKS 1510-089, reported in Fig. 124"! show that during 
the flaring state, the IC peak dominated over the synchrotron one by more than one order of 
magnitude. Indeed, even if both the synchrotron and IC peak frequencies are not sampled in our 
data set, we can estimate their typical peak flux by extrapolating the LP best fit model for the 
7-ray data, and by fitting the radio, optical, and UV SED points by means of cubic function. In the 
case of the synchrotron component, simultaneous mm and optical/UV data indicate that the peak 
flux should not be higher than a few times 10 _11 erg cm -2 s _1 . In the case of the IC component, 
since the spectrum is mildly curved with the peak energy below the threshold of our analysis (200 
MeV), the peak flux must be a few times 10 -10 erg cm~ 2 s _1 . The presence of a pronounced BBB 
and the Compton dominance of about 10, during the flaring states, suggest that the contribution 
from photons originating outside the jet is appropriate to model the broadband SED. 

The absence of X-ray/7— ray nux correlation hints that the ERC flux variations depends on 
a change in the high energy spectral index of the electron distribution, instead of a change in the 
external radiation field. In Sect. 14.21 we noted that the lack of 7— ray/UV correlation suggests that 
BBB variations are excluded as a main driver of the 7— ray variability. We tested the change of 
electron distribution by studying the correlation between the logarithms of the 7— ray and of the 
Optical R fluxes (see Fig. [T6|) . The relation between the 7-ray and optical flux, fitted by means of 
power law, returns an exponent of ~ 0.5 (or ~ 0.6 if we use the 13-day shifted optical light curve 
discussed in Sect. 14. 2p . Despite the large scatter, this value is very interesting, because synchrotron 
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and IC fluxes correlate very differently in the case of SSC, ERC/BLR, and ERC/DT. Moreover, the 
correlation dep ends on the low and high energy range chosen respectively for the synchrotron and 



IC component (jKatarzyhski et al.ll2005l ). To have an estimate of the power-law exponent for this 



correlation, we reproduced nu merically a set of SEDs comparable to that observed for PKS 1510 



089, using a well tested code (jTramacerel 120071 ; iTramacere et al.l J2009I : iTramacere fc Tostil 12003 ) 



We calculated an ERC-dominated model, where the 7— ray emission is dominated by the BLR 
contribution and the electron distribution is a BPL. We fixed all the parameters and we increased 
the high energy index of the electron distribution from -3.5 to -2.5 to evaluate the correlation 
between the optical energy emission (uF(u) at 10 14 Hz) and the integrated F(E > 200 MeV) 7— ray 
flux, taking separately into account the three IC components SSC, ERC/BLR and ERC/DT, as 
shown in Fig. [25j The resulting flux relations were fitted by simple power laws and the resulting 
exponents were 3.1 for the SSC, 1.8 for ERC/DT and 0.4 for ERC/BLR. The last one was the only 
one found to be consistent with that observed for PKS 1510-089 (~ 0.5-0.6) and therefore, we can 
disfavor both SSC and ERC/DT. 

The SEDs of the ERC/BLR model in the upper panel of Fig. [25] have a curvature more 
pronounced than the ERC/DT ones. This is due to the Klein-Nishina cross section, because UV 
photons emitted from the accretion disk and reflected towards the jet by the BLR are blue shifted 
roughly by a factor of V. For V ~ 10, the typical energy of these photons in the emitting region rest 
frame is then ~ 10 16 Hz, hence the IC scattering with the electrons with 7 ~ 1000 occurs under a 
mild KN regime. The resulting smooth curvature has a value of ~ 0.1 that is compatible with the 
observed one. 



5.2. Spectral energy distribution modelling and jet energetics 

We attempt a leptonic ERC/BLR-oriented SED modelling. We aim to reproduce the SEDs 
for the 3 flares with simultaneous data from radio to 7— ray energies, namely flares b,c, and d. 
Moreover we try also to fit the quiescent state. Although variability timescales reached values of 
a fraction of day, LAT data required longer integration times to produce a SED enough good for 
spectral modelling. The intermediate fluxes of the synchrotron and BBB components, observed 
during the 7— ray integration period (red filled circles in Fig. [2l"|) . are used in the fitting procedure 
as representative of the flare- averaged state of the low energy bump. 

Since the correlation between the 7— ray and the optical fluxes (see. Sect. 15.11 and Fig. I25j) 
favours an ERC/BLR scenario, we assume that the dissipation zone is in the sub-pc scale. This 
is also consistent with the mild curvature observed in the 7— ray spectra. Indeed, as shown in 
the previous Section, the ERC/BLR process occurs under the KN regime leading to the curved 
MeV/GeV spectral shape that matches the one observed in the Fermi spectra. 

We assume a jet viewing angle of 6 = 2.5° for both the flaring and the quiescent states. During 
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the flaring states, we choose a bulk Lorentz factor in the range [14-16], resulting in a beaming factor 
range of [20-21.5] that is compatible with the VLBI observations. During the quiescent state, we 
use a bulk Lorentz factor of 12.0, corresponding to a beaming factor of about 18.8. 

As a further step, we estimate the accretion disk physical characteristics using the UV data. 
We use the UV observations during the lowest synchrotron state of our dataset as an estimate for 
the upper limit of the accretion disk luminosity (X<f). According to the observed UV flux and to 
the luminosity distance, we set a reference y alue of L^ ~ 3 x 10 45 erg/s. We mod el the accretion 
disk, following the prescription in lKineJ ((2008) and rearranging the expression as in lGhisellini et al 
(|2009l ). using a multi-temperature blackbody with a temperature profile given by: 



disk 



(R) 



ZR S L 



1 



(m n 



(2) 



167re<7sB.R 3 

where osb is the StefanBoltzmann constant, R$ = 2GMbh/c 2 is the Schwarzschild radius for a 
black hole (BH) mass Mbh, 3Rs is the last stable orbit in the case of Schwarzschild BH, and e is the 
accretion efficiency that is linked to the bolometric luminosity and to the accretion mass rate M, 
by L<i = eMc 2 . We assume that the radiative region of the disk extends from ~ 3Rs to ~ 500Rs ■ 
Since Tdisk(R) peaks at R ~ ^Rs, we can use the UV data to constrain the peak of Tdiks(R) and 
we can use Eq. [2] to estimate Rs- From our data we get a value of Tj 



peak 



disk 



4 x 10 K. Assuming 
as a reference value for the accretion efficiency e ~ 0.1 we get Rs — 1.6 x 10 14 cm. The BH mass 
Mbh is then about 5.4 x 1O 8 M and the accretion rate of about 0.5 M /year that corresponds to 
~ 0.04 times the Eddington accretion rate f M^). The va l ue of the BH mass estimated from our 
UV data is compatible with that obtained bv lOshlack et al.l (|2002j ) using the virial assumption with 
measurements of the H/3 FWHM and luminosity (M BH ~ 3.9xlO 8 M ), and by Ixie et aD (J20od ) 
(Mbh =* 2xlO 8 M ). 



The radius of the BLR can be estimated using the empirical relation by iKaspi et al.1 ((2005) 
and bentz et all (J20od ) : 



Rblr =* 20 
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.10 44 erg s 

where v* ~ 5.878 x 10 14 Hz. According to the UVOT data we estimate v*L(v*) < 2 x 10 45 erg/s, 
and we get a value of Rrt.fi < 2.3 x 10 cm . We use BLR radius of about 1.6 x 10 17 cm with 
reflectivity value of 0.1 (jGhisellini et al.ll2009r). The DT is model led as a BB with a temperature 
Tb,t of about 100 K, and a reflectivity of 0.3 (jGhisellini et al.ll2009l ). The distance is set to a typical 
value of ~ 10 18 cm, and it is fine tuned in the fit in order that the BLR/DT+BLR/ERC correctly 
match the X-ray and hard X-ray data. 

We use an emitting region size (R ra d) i n the range [2 — 5] x 10 16 cm, that is almost twice 
the value estimated from the fast variability. In this regard, we note that since we are describing 
flare-averaged states, with integration times of the order of a few weeks, the discrepancy with the 
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fast variability estimate is not problematic. As electron energy distribution N(j) we use a broken 
power-law: 

iV( 7 )oc( 7 " fOT ^n<l<^ (4) 

[j pl for 7fer < 7 < Imax 

with the number density of the electrons: 

N(j)dj. (5) 

,min 

The low energy spectral index (p) is chosen ~ 1.9, as hinted by the typical X-ray photon index (see 
Sect. 15. ip : the value of the break energy (jbr) is chosen to match the position of the peak energy 
of the IC bump. It is of the order of [250-300] and it is tuned in the three flares to fit the data. 
The high energy spectral index (pi) is chosen according to the average spectral index of the LAT 
data, and the reference value is ~ 3. Magnetic field intensity B is chosen to be 1.0 G for flares b 
and c, while for the flare d we need to use a larger value of B and a more compact region size. 

All the flares are assumed to occur at dissipation distance from the disk (Rdi. ss ) that is within 
the BLR and the DT (Rdiss < Rblr, Rdiss < Rdt) ■ The resulting best-fit of the MW SEDs are 
reported in Fig. [2^| and the corresponding values of the best fit parameters for the three flares are 
reported in Tab. [5j 

To have an indication of the change in the energetic contents of the jet as a function of the 
flaring activity, we try to fit also the quiescent 7— ray SED. Considering the lack of MW data during 
this time interval, we assume that the X-ray SED is close to the state observed during flare b, as 
supported by the low X-ray variability. From the MW light curves reported in Fig. we infer 
that the optical/UV flux is at a level of few times 10 -12 erg cm~ 2 s . We can have an acceptable 
description of the MW SED during the quiescent 7— ray period using parameters comparable to 
those chosen for the flare b, but we need to decrease significantly 7^ to ~ 65. The best-fit model 
is represented by the black dash-dotted line in the top panel of Fig. 



It is interesting to study the evolution of the energetic content of the jet resulting from the 
SED modelling discussed above. We evaluate the total kinetic power of the jet as: 

L jet = 7rR 2 rad cT 2 (U' e + U' B + U' p ) (6) 

where the rest frame magnetic energy density is given by U' B = B 2 /8tt, the electron energy density is 
C/Jj, and Up = Q.l N P m p (? is the cold -proton energy density, assuming that we have one cold proton 



per ten electrons (jSikora et al.ll2009l ). The power carried by the jet in terms of radiation is given by 



W ' = -CT 7m eC 2 7V( 7 )d 7 
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L ra( i ~ L'r 2 /4. We evaluate L' summing up the numerical integration of each radiative component 
(synchrotron, SSC, ERC/BLR, ERC/DT) as observed in the jet rest frame. In Tab. [6j we report 
the results for flare b, c, d and the quiescent state, corresponding to the physical parameters reported 
in Tab. [5j 

The total kinetic power of the jets is almost steady, except during the flare d when it increased 
by a factor of 2. Indeed, for the case of flare b and c, we find values of Lj e t that are comparable 
with that of the quiescent state. On the contrary, we note that the value of the electron energy 
density during all the flaring states is much larger than that estimated during the quiescent state, 
and, during the flare d, it increased with respect to the quiescent state by about a factor of 2.7. 
The low values of U' e /U' B reflect the low SSC contribution. Indeed, for the choice of our model 
parameters, the SSC is always negligible compared to the other radiative components. As final 
remark, we note that the radiative efficiency of the jet (L ra d/Lj e t), is ~ 0.03 during the quiescent 
state, and increases up to ~ 0.2 during flare c. 

We now compare the jet total kinetic l umino sity with the accretion disk luminosity. According 



to the analysis reported in lGhisellini et al.l (|2010l ). if the jet power comes from the accretion process, 
then the accreting mass has to account both for the disk luminosity and for jet power. We can 
write: 



etotMc 2 = e D Mc z + e jet Mc 2 . (7) 

The ratio of Lj et to Lp, reported in Tab. [6]is equal to tjet/tD- Using as eo the typical value of 0.1 
we need a total efficiency etot ^ 0.26, that is still compatible with the maximum efficiency in the 
case of a Kerr BH. 



5.3. Comparison with other Fermi-hAT FSRQs. 

It is interesting to compare the flaring activity of PKS 1510-089 with other FSRQs observed by 
Fermi. 3C 454.3 (z=0.859), reached an isotropic lumi nosity above 100 Me V L(E > 100) ~ 8 x 10 48 



erg/s during the period of 2008 August/September (lAbdo et alJl2009al ). It became the highest 



ever recorded object in the gamma rays during the flare in 2009 December (JEscande fc Tanaka 



2009), with a luminosity L{E > 100) ~ 1 x 10 49 erg/s, about five times larger than the maximum 
flare resolved luminosity of PKS 1510-089. The source during the flares exhibited rapid variability 
down to sub-daily timescales, that is comparable to that observ ed in PKS 1510-089 . The 7— ray 



spectra of 3C 454.3 object show a spectral break around 2 GeV (jAbdo et alj|2009al ). In this case 
we see a relevant difference w.r.t. to PKS 1510-089, that exhibited a mild curvature spectrum. A 
possible reason, as explained in the previous Section, is that in the case of 3C 454.3 the dominant 
external photon field originates in the DT and not in the BBB. Thus the IC process occurs under 
the TH regime, and the break reflects the break in the electron distribution. The black hole mass 
is estimated to be Mbh — 4 x 10 9 M q , and the disk luminosity of the order of Ln — 2 x 10 46 
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erg/s, both these values are one order of magnitude larger than PKS 1510-089. Another interesting 
difference is that the 3C 454.3 exhibits no BBB features during the high optical and 7— ray outburst. 



The FSRQ PKS 1502+106 (2=1.839), during the outburst in 2008 August (jAbdo et alJl2010d ). 
reached a luminosity L ~ 1 x 10 49 erg/s, showing a very large Compton dominance up to 100. As in 
the case of 3C 454.3, this object exhibited no BBB feature at UV energies. The flare of this object 
was an isolated episode, rather than a flaring sequence. The 7— ray spectrum showed a curved 
shape, similar to PKS 1510-089, during the full integration period, and during the post-flare state. 
The curvature values, 0.1-0.2, are comparable to those observed in PKS 1510-089, hence also in 
this case an origin in the KN effect is possible, as well as a curved electron distribution with the 
IC process in TH regime. The MW variability of this object is very different from that observed 
in PKS 1510-089. Indeed, the X-ray flux followed both the 7— ray and the optical flux. Moreover, 
there was a strong correlation between the UV and the 7— ray. Probably this can be explained with 
a sharp change in the density of the radiating electrons, or in the beaming factor. 

PKS 1454-354 (2=1.424), another di stant FSRQs discov ered in the 7— ray by Fermi during the 



2008 August/September flaring activity (JAbdo et al.ll2009bl ). reached a luminosity L(E > 100) ~ 



7 x 10 48 erg/s, and exhibited a rapid variability with hours/days timescales. 

We conclude that the most rapid timescale for the 7— ray variability (down to 6/12 hr) is 
common to this class of objects, despite the bolometric 7— ray luminosity. This timescale is limited 
by the minimum LAT integration time required to extract a flux for object with F{E > 200MeV) 
of the order of 10 -5 ph cm~ 2 s _1 , hence a faster variability could be possible. Compared to those 
listed above, PKS 1510-089 seems to be a less powerful object in terms of both 7— ray luminosity 
and thermal luminosity. The curvature in the 7— ray spectrum, when compared to the spectral 
break observed in other sources such as 3C 454.3, may be understood as a signature of the KN 
regime, hence of a dissipation region within the BLR. 



6. Conclusions 

We presented MW observations of PKS 1510-089 during a period of about 11 months, when the 
source exhibited a strong evolution of its broad-band SED, characterized by a complex variability 
both at Optical/UV and 7— ray energies, with timescales detected down to the level of 6/12 hr. 
The 7— ray flux shows the usual harder when brighter trend, for flux levels of F(E > 0.2GeV) > 
2.4 x 10 -7 ph cm~ 2 s _1 . For lower flux levels the trend seems to be the opposite, but the low statistics 
during these dimmer states does not allow us to give a purely physical interpretation. The spectrum 
shows a mild curvature, both during the flares and during the full-period integration, well described 
by a log-parabolic law, and can be understood as a signature of the KN effect. 

The 7— ray flux shows a complex correlation with the other wavelengths. There is no correlation 
at all with the X-ray band, a week correlation with the UV, and a strong correlation with the optical 
(R) flux. Thanks to the unprecedented continuous LAT 7— ray sky survey, we were able to find a 
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lag in the 7— ray light curve, with the 7— ray band leading the i?-band by about 13 days. Because of 
this complex multi-band variability, we assume that a change in the beaming factor cannot account 
for the flaring activity of this object. Indeed, assuming as the main driver a change in the beaming 
factor, we would expect at least a weak correlation between the X-ray and the optical band, or 
between the X-ray and the 7— ray. This is also in agreement with the absence of correlation between 
the jet kinematic power(Lj e t) and the flaring episodes observed in our data set. 

UV data allowed us to estimate the mass of the BH of ~ 5.6 x 1O 8 M0. This value, that is 
in agreement with other estimates based on different methods, is smaller compared to the very 
luminous and distant FSRQs, with BH masses of the order of 10 9 Mq. As a consequence of the 
estimated BH mass and thermal component luminosity, also the accretion rate of ~ 0.04Me^, is 
lower when compared to the expectation in case of FSRQs (~ O.lMEdd)- 

Due to the low redshift of the source, the bolometric isotropic 7— ray luminosity is also smaller 
compared to other distant FSRQs observed by Fermi. Indeed, PKS 1510-089 has a typical 7— ray 
luminosity and BH mass about one order of magnitude lower compared to sources like 3C 454.3 or 
PKS 1502-106. 

Despite the relatively lower power in the thermal and non-thermal outputs, PKS 1510-089 ex- 
hibits a quite large Compton dominance, as observed in the most powerful FSRQs, and prominent 
a BBB signature. The object could be a representative of an aged FSRQ, hence the analysis here 
presented is relevant in order to understand the evolution of these objects. 

We note the puzzling feature of the BBB UV shape. The BBB was still prominent during the 
historical maximum optical state during 2009 May, although the optical/UV spectral index was 
softer compared to that in quiescent state. 

The analysis presented here shows the importance of the MW monitoring of blazars, inde- 
pendent of the 7— ray triggering. Indeed, only comparing the flaring and quiescent states, and 
understanding the evolution of the physical parameters as a function of the flaring activity, is it 
possible to discriminate among the possible physical scenarios. 
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Table 1: Flaring activity of PKS 1510-089 from August 2008 until June 2009 



source state date start date stop MJD start MJD stop 



flare a 


2008-08-30 


2008-09-26 


54708 


54735 


quiescent 


2009-09-30 


2010-01-01 


54739 


54832 


flare b 


2009-01-04 


2009-01-27 


54835 


54858 


flare c 


2009-03-10 


2009-04-09 


54900 


54930 


flare d 


2009-04-15 


2009-05-12 


54936 


54963 
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Table 2: Unbinned likelihood LRT summary 



time range 


-LogLikc(PL) 


-LogLikc(LP) 


-LogLikc(BPL) 


LRT(PL/LP)/P.(LP)** 


LRT(LP/BPL)/P.(BPL)** 


full 


323082.6 


323056.4 


323062.8 
323059.0* 


56.2/< 10~ 5 % 


-12.8/NULL 
-5.2/NULL 


quiescent 


83057.7 


83057.5 


83057.6 


0.4/47.3% 


-0.2/NULL 


flare a 


28908.9 


28902.4 


28903.5 
28902.5* 


6.5/99.97% 


-2.2/NULL 
-2.0/NULL 


flare b 


23932.7 


23928.5 


23927.9 
23927.9* 


8.4/99.62% 


1.2/72.67% 
1.2/72.67% 


flare c 


38328.9 


38318.5 


38320.2 
39319.4* 


20.8/99.99% 


-3.4/NULL 
-1.8/NULL 


flare d 


31326.1 


31322.0 


31323.3 

31322.8* 


8.2/99.58% 


-2.6/NULL 
-1.6/NULL 



* BPL fit by means of loglikelihood profile 

** P.(LP) and P.(BPL) are the cumulative distribution functions of the LRT statistics, evaluated at the 
LRT value actually observed. These probabilities are evaluated using as reference distribution a \d 
distribution with the number of degree of freedom (d) equal to the difference in the number of free 
parameters in the two models. 



Table 3: Unbinned likelihood spectral fit results: 



time range 


Oiy 


PL 
F 100 


-LogLikc 


a-^ 


LP 


F inn 


-LogLikc 


a-y 


a -n 


BPL*** 


F inn 


-LogLikc 








full 


2.44±0.01 
2.43±0.07 
2.55±0.01 
2.35±0.04 
2.37±0.3 
2.44 ±0.04 


1.32±0.03 
0.27±0.03 
1.55±0.02 
2.1 ±0.1 
3.9 ±0.1 
2.9 ±0.1 


323082.6 
83057.7 
28908.9 
23932.7 
38328.9 
31326.1 


2.23±0.02 
2.3±0.1 
2.2±0.1 
2.0±0.1 

2.13±0.06 
2.24±0.08 


0.09 ± 0.01 
0.03 ± 0.05 
0.19 ± 0.06 
0.10 ± 0.03 
0.1 ± 0.02 
0.09 ± 0.03 


1.12±0.03 

0.25±0.04 

1.16±0.12 

1.7±0.1 

3.3±0.2 

2.5±0.1 


323056.4 

83057.5 
28902.4 
23928.5 
38318.5 
31322.0 


2.30±0.02 

2.5±1.5 
2.39±0.01 
2.26±0.05 
2.28±0.04 
2.31±0.08 


2.63 ± 0.04 
2.4 ± 1.5 

3.22 ± 0.04 
3. ±0.03 
2.9±0.2 
2.6 ±0.1 


980±130 


1.37±01 


83057.6 


quiescent 


700±8000 


0.3±02 


28903.5 


flare a 


1500±40 


1.37±0.01 


28903.5 


flare b 


3400 ±800 


1.9±0.1 


23927.9 


flare c 


1918 ±593 


3.7 ±0.1 


38319.9 


flare d 


1000 ±300 


2.6 ±0.1 


31323.3 



10 _q ph cm - ' s~ J 
* McV 

BLP fit by means of loglikclihood profile 



CO 
00 
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Table 4: Spectral analysis of XRT data. In the last column we report the reduced \ 2 an( i, in 
parentheses, the degree of freedom. 



Observation Date 



Start Time 
MJD (days) 



Norm 
10- 4 



a x 



flux 0.3-10 keV 
10~ 12 erg cmr' 1 s~ 



Xl/doi 



2009-01-10 
2009-01-11 
2009-01-13 
2009-01-14 
2009-01-16 
2009-01-16 
2009-01-21 
2009-01-25 
2009-03-06 
2009-03-11 
2009-03-12 
2009-03-17 
2009-03-18 
2009-03-19 
2009-03-20 
2009-03-22 
2009-03-22 
2009-03-23 
2009-03-24 
2009-03-25 
2009-03-26 
2009-03-27 
2009-03-28 
2009-03-30 
2009-04-04 
2009-04-10 
2009-04-27 
2009-04-28 
2009-04-29 
2009-05-01 
2009-05-02 
2009-05-03 
2009-05-04 
2009-05-05 
2009-05-07 
2009-05-11 
2009-05-12 
2009-05-13 
2009-05-14 
2009-05-30 
2009-06-07 
2009-06-13 
2009-06-20 



54841.8 


9 4+O6 


1 or+0.07 

1 - oa -0.08 


8 5+°- 8 


54842.9 


q 4+0-6 


1 ao+0.07 
1 - zo -0.07 


9 4 +0 -I 


54844.0 


9.9+83 


1 46+ - 08 


7 q+o-e 


54845.8 


10-4±8;S 


1 41 +0.07 
i - 41 -0.07 


O+0.7 
0*0-0.9 


54847.7 


10.2+83 


1 AK+ om 
i - 4a -0.07 


0+O.8 
°- z -0.8 


54849.2 


8.o±8:S 


1 44+ - 12 

1 - 44: -0.12 


6 5+°- 8 

D -°-0.8 


54852.8 


0+0.6 

y -°-0.6 


1 41+007 
i - 41 -0.07 


O+0.7 
».d_0.7 


54856.8 


9.6I°; 7 r 


1 S6+ 009 

l.OU_ 09 


7+O.8 
°*'-0.8 


54896.9 


8.9±M 


1 r+0.2 
i,o -0.2 


7 1+ 1 ' 2 

'•-■--1.4 


54901.6 


ii-s±g:S 


1 62+ 006 

1 - uz -0.06 


7 0+O.6 
'•O-0.5 


54902.6 


13 1+ 06 


1 ^8+ - 05 

1 -°°-0.06 


q 1+0-6 


54907.2 


10 7 +0 - 6 


1 fi +0.07 
J "°-0.07 


7 o+0-6 
' - z -0.6 


54908.0 


O q+0.5 
y - y -0.5 


1 AK+0-06 
i - 4a -0.06 


0+O.6 
O-0.5 


54909.8 


12 l +0 - 8 


1 RQ+0.09 
l.U»_ .09 


7-418:1 


54910.9 


8.711 


1 O+0.14 

i -°-0.14 


8.511J 


54912.1 


q q+0-9 


1-4618:11 


7.918; 8 


54912.1 


16.5^1 


i-69l8:l 


10.1113 


54913.5 


9-s±g:? 


i.47l8:l 


7 0+0.6 

'•°-0.7 


54914.1 


10 2 +0.9 


1 5fi+ - 12 


7 0+O.8 
'•°-0.7 


54915.6 


11 1+ ' 8 


1 ci +0.08 
i - ol -0.08 


84 +0.9 


54916.7 


10.6+} 


1 CJ4+0.12 

r-o / i_o.i2 


7 6+1-1 

'•"-1.3 


54917.2 


9.9+1 


1 ro + 0.12 

1 - oo -0.12 


70 + I.I 
'•°-0.8 


54918.2 


10.0I8:? 


1 or+0.08 

1 - oo -0.08 


q 1+0.9 

y * ± -0.7 


54920.4 


11.818-J 


1 41+0-08 
i -* i -0.08 


10+I 


54925.5 


8.H13 


1 0+0.2 
i, ^ 1 -0.2 


81 2 


54931.0 


10.518-1 


1 e+0.07 


0+0.6 

O-0.7 


54948.6 


12.711-j 


1 62+ 013 

1 - DZ -0.13 


8-411.3 


54949.6 


12.0+2 


1 -1O+0.13 
1 - lo -0.13 


15+i 


54950.8 


11.711+ 


1 36+ 011 

1 - oo -0.11 


10.5113 


54952.6 


IO.OI2 


1-5183 


77 +i.7 
'•'-1.2 


54953.4 


13.0+1 


1-5818:1 


q+1.2 

y -0.8 


54954.8 


11.911-j 


1 tc+0.2 

i *°-0.2 


913 


54955.7 


10.0+2 


1 0+0.2 
s — c -0.2 


"12 


54956.7 


io.5i8:8 


1 - oo -0.12 


7 7+0.7 
'•'-0.8 


54958.4 


io.ii8:l 


1 90+0. 08 
1 - zo -0.08 


io.i+;- 9 


54962.5 


5-4111 


q+o- 3 

u - y -0.3 


io±i 


54963.3 


13.6+83 


1 4^+0-06 


11 2+ - 8 

iJ -' z -0.8 


54964.1 


19 q + 0.6 

iz -°-0.6 


1 £+0.06 
i - u -0.06 


O+0.5 

o.o_ 06 


54965.3 


11 1+ - 6 

ii,1 -0.6 


1 AS+0-06 
i - 48 -0.06 


8.618:1 


54981.4 


12.011 


1 7+0-4 
L' -0.5 


7+3 
'-2 


54989.7 


17.0+^ 


1 8+0- 2 
•■■■O-0.2 


9.2113 


54995.6 


".611:1 


1 or+0.14 
1 .°°-0.14 


10.6113 


55002.5 


16.0113 


i.86l8:l 


0+0.8 
o-^-o.s 



0.697 


43) 


0.817 


46) 


0.978 


38) 


0.811 


46) 


0.778 


42) 


0.571 


20) 


0.898 


51) 


0.832 


30) 


0.980 


11) 


1.067 


68) 


0.897 


78) 


1.184 


61) 


0.949 


59) 


1.049 


35) 


1.210 


20) 


1.325 


30) 


0.506 


20) 


0.983 


30) 


0.770 


23) 


1.293 


37) 


1.084 


19) 


0.840 


19) 


0.777 


32) 


0.855 


31) 


1.775 


12) 


1.098 


46) 


0.502 


15) 


0.446 


13) 


0.951 


25) 


1.327 


08) 


1.241 


29) 


0.838 


10) 


1.428 


08) 


0.838 


22) 


1.122 


30) 


1.660 


05) 


0.931 


62) 


0.822 


59) 


1.175 


56) 


0.768 


02) 


1.575 


16) 


1.185 


20) 


0.888 


24) 
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Table 5: Best-fit parameters for the SED modelling in the case of BPL electron distribution. 



flare L D M B h Rblr Rdt T R ra d B N e ^ min ^ max -y br p pi 

10 45 crg/s 10 8 M e 10 17 cm 10 17 cm 10 17 cm G #/cm 3 

~~b 3^0 5J L6 12~F 

c 



quiescent 



.0 


0.32 


1.0 


550 


1.0 


2.2xl0 4 


220 


1.95 


3.15 


.5 


0.32 


1.1 


300 


1.0 


2.2xl0 4 


280 


1.80 


3.15 


.0 


0.25 


2.2 


800 


1.0 


7.0 xlO 3 


200 


1.90 


3.20 


.0 


0.43 


1.0 


350 


1.0 


2.2xl0 4 


65 


1.95 


3.20 



Rdiss<RBLR, Rdiss<R D T, T^ s f=4xl0 4 K, T DT = 100K, 9 = 2.5 C 



Table 6: Jet power for the case of BPL modellig. 



flare 


Lj e t 

10 45 crg/s 


l^rad 

10 45 crg/s 


Lei Ljet 


Lb /Lj e t 


Lrad/ Ljet 


U e 

crg/cm' 


U'elU' B 


Ljet/Ln 


b 


2.32 


0.26 


0.025 


0.316 


0.112 


0.0031 


0.08 


0.77 


c 


2.17 


0.46 


0.027 


0.502 


0.211 


0.0026 


0.05 


0.72 


d 


4.83 


0.34 


0.016 


0.606 


0.078 


0.0049 


0.03 


1.61 


quiescent 


2.31 


0.07 


0.016 


0.423 


0.031 


0.0015 


0.04 


0.77 



a 3 

o 

o 



2 - 



4x10 - 



3x10 



> 

o 

s 

o 

o 

(N 
A 

S 

3 



2x10 



1x10 
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680 700 720 740 760 780 800 820 840 860 880 900 920 940 960 980 1000 1020 

Time (MJD - 54000) 



Fig. 1. — Upper Panel: the 7— ray photon index (a 7 ), as a function of the time, for weekly and 
daily binning. In case of daily binning only observations with a Test Statistic > 10 are taken into 
account. The Test Statistic is defined as TS: = -2 Log(L0/-Ll), where LI and LO are the likelihood 
when the source is included or not. Lower Panel: weekly, and daily( TS>10) fluxes light curve. 



2x10" 



> 



o- 

O t« 

A 'a 1x10 

m 5 



~r 

♦ — ♦ GASP r *(rescaled) 



37 



i i I i i i i i I i i i i i I i i i i i I i i i i i 




3x10 



705 710 715 720 725 730 735 

Time (MJD - 54000) 



740 



> 

©'</> 2x10 

OcM 

8^ -6 

^3 1x10 







♦ — ♦ GASP r *(rescaled) 



i i I i i i i i I i i i i i I i i i i i I i i i i i 




830 835 840 845 850 855 860 865 

Time (MJD-54000) 



5x10 
% 4xl0" 6 



s^r 3xio"° 

^ 'd 

W ^ 2x10 

x 3 

E lxio" 6 



■ 



♦ ♦ GASP r *(rescaled) 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I i- 




890 895 900 905 910 915 920 925 930 935 940 

Time (MJD-54000) 



£■„ 4x10" 



i i i i i i i i i i i i i i i 



i I i i i i i I i i i i i I i i i i i I i i i i i -I 




940 945 950 955 960 965 970 

Time (MJD-54000) 



Fig. 2. — from top to bottom, flare a, b, c, and d showing 1-day binning (red points) and 12 hr 
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Fig. 3.— LAT SED of PKS 1510-089 extracted for the full period (black points) and for the 
quiescent state (green points). The red upward arrow indicates the highest energy event within 
95% of the PSF for the whole period data set. The dotted line, the dashed line, and the dot- 
dashed line, represent the best fit model of the full period by means of PL, LP, a BPL distribution 
respectively, with uncertainties. The residuals in the lower panel refers to the PL model. 
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Fig. 4. — LAT SED of PKS 1510-089 extracted for the flaring states and for the quiescent state 
(green points). 
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Fig. 5. — Upper Panel: weekly and daily scatter plot of the Flux( E> 200 MeV) vs. the photon 
index (TS > 10). The inset shows the same for Flux( E> 200 MeV) > 2x10 ph cm~ 2 s _1 . Lower 
Panel: Scatter plot of the number of photons predicted by the best-fit model vs. the photon index, 
for weekly and daily integration (TS >10). 
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Fig. 6. — Upper Panel: The histogram of the photon index for, a weekly integration time, before 
(blue) and after (red) MJD 54905, respectively. Upper Panel: The same as in the upper panel, in 
the case of daily integration. 
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Fig. 7. — The MW light curves, from 2008 April to 2009 June. The vertical dashed lines show the 
the four flaring episodes and the quiescent state. 
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Fig. 8. — Scatter plot of the flux in 0.3 — 10.0 keV range vs. the X-ray photon index (ax)- The 
dashed lines represent a linear bet fit model. Using the Monte Carlo method described in Sect. 
12.2. H we get a correlation coefficient r = —0.31 with a 95% confidence limit of —0.55 < r < —0.05 
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Fig. 9. — XRT SEDs averaged during the flare b, c, d and duirng the post d flare period. The 
violet boxes represent the BAT spectrum averaged over 5 years 
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Fig. 10.— Scatter plot of the XRT flux in the 0.3-10 keV band vs. the Fermi-LAT flux (E> 200 
MeV). 
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Fig. 11. — Radio to UV SEDs built using only data simultaneous within a daily time span. 






> 



0.15 



0.1 



0.05 



M 

o 



-0.05 - 



-0.1 



-11.5 -11 

log(vF(v)) R filter 



1 






1 






- 


| • UVW2/B| 




• 


• 










- 






• • 






- 


- 






• 
• 


• 
• 

• 




- 


- 








• 


• 


- 


- 




• 
• 






• 


— 


- 








• • 




- 


- 








• 


• 


- 


- 










• 


- 


1 






1 






- 



-10.5 
(erg cm s ) 



Fig. 12.— The ratio of uY(v) in the UVOT UVW2 filter to ^F(z^) in the UVOT B filter, as a 
function of vF(v) in the optical R filter. 
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Fig. 13. — Panels a,b: UVOT light curves. Panel c: Hardness ratio of the UVOT spectra 
(UVW2/UVF) evaluated as a function of the time. 
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Fig. 14.— Scatter plot of the Fermi-LAT flux (E> 200 MeV) vs. the UVOT uF{v) in the UVW2 
filter. The red dashed line represents the best-fit power law model. The correlation coefficient of the 
logarithm of the UV and 7-ray fluxes is r = 0.2 with a 95% confidence interval of 0.05 < r < 0.34. 
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Fig. 15. — Optical and near-IR light curves from the GASP project. 
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Fig. 16.— Scatter plot of the Fermi-LAT flux(E> 200 MeV) vs the vF{u) in the optical R band. 
The dashed line represents the best fit power-law model. The correlation coefficient of the logarithm 
of the optical and 7-ray fluxes, evaluated through the Monte Carlo method described in Sect. I2.2.1[ 
is r = 0.42 with a 95% confidence interval 0.36 < r < 0.46. There is an evident dispersion in the 
scatter plot, probably due to the inter-band time lags showed in Sect. 14.21 
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Fig. 17. — Top panel, left: DCCF for the whole analysis period, the shaded box concerns the 13 
days lag. Top panel, right: fit by means of Gaussian distribution of the 13 day lag shaded in the 
left panel, the fit returns a lag of -13. 4± 0.2, gamma leading optical. Middle panels: the DCCF 
for MDJ<54890 (left panel) and MDJ>54890 (right panel). The lag of about 13 days with the 
7-ray leading the optical band is still there. Bottom panels: a lag of about 13 days (gamma leading 
optical) is also present for flare d and b individually. 
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Fig. 18. — Upper Panel: the optical (R) light curve, obtained shifting the time according to the 13 
days 7— ray lag reported in Sect. 14.21 Lower Panel: The Fermi-LAT light curve 
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Fig. 19.— Scatter plot of the Fermi-LAT fiux(E> 200 MeV) vs the vF{y) in the optical R band, 
using the time-shifted fluxes reported in Fig [18] . The dashed line represent the best fit by means 
of power-law model. The correlation coefficient in this case increases to r = 0.62 with a 95% 
confidence interval 0.57 < r < 0.66. 
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Fig. 20. — Upper panel: The 15 GHz MOJAVE VLBA parsec-scale core light curve. Lower panel: 
GASP radio light curve. 
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Fig. 21.— Stokes I CLEAN images of PKS 1510-089 observed by VLBA on 2009 April 9, at 15, 
24, and 43 GHz. The lowest contour and peak intensity are 0.7 mJy/beam and 1.47 Jy/beam 
(15 GHz), 0.7 mJy/beam and 1.56 Jy/beam (24 GHz), 2 mJy/beam and 1.85 Jy/beam (43 GHz). 
Contours are plotted with a step x4. Natural weighting of visibility data is used, HPBW beam 
size is shown in the lower left corner. Angular size of 1 mas corresponds to 5 pc. 
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Fig. 22. — Spectral index s r (F(v) oc v Sr ) map between 15.4 and 23.8 GHz (shown in color) of 
PKS 1510-089 as observed by the VLBA on 2009 April 9. The overlaid contours represent total 
intensity at 15.4 GHz (see Figure [21] for details). The spectral index map was smoothed by a 
median filter with a 0.6 mas radius. 
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Fig. 23. — Optical R — H spectral index vs. the LAT Flux above 200 MeV. The black points refer 
to the 2008 August to 2009 February period, and are poorly sampled to be used to investigate the 
optical/7— ray connection. The blue points overlap mainly the flare c, in particular the subflare 
peaking at about MJD 54909 (see fig. [2j panel d). The red arrows shows the chronological sequence. 
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Fig. 24. — Red solid circles correspond to simultaneous Optical/UV data representing the in- 
termediate state during the 7— ray integration period. Blue and green circles correspond to the 
highest and lowest flux state, respectively, observed during the same integration period. The 
Swift-XRT data (orange points) are integrated during the same interval of the LAT data, and the 
Swift-BAT data represent the 5-year fluxes discussed in Sect. 14.11 . The black squares represent 
GASP radio data, integrated during the flares. The cyan points correspond to the LAT quiescent 
state SED. The magenta points represent the 7— ray SED integrated during the flare. The blue, 
green, and red dotted lines represent the PL best fit of the 7— ray data, for a daily integration and 
simultaneous to the Optical/UV data with the same color. The thin black dotted line represents 
the flaring synchrotron emission. The thin black dashed line corresponds to flaring SSC emission. 
The red and blue thick dotted lines correspond to the dusty torus and BBB emission respectively. 
The red and blue thin dashed lines correspond to the ERC/DT and ERC/BLR flaring emission 
respectively. The solid black thick line represents the sum of all the flaring model components. The 
thick dot-dashed line represents the sum of all the model components for the quiescent state alone. 
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Fig. 25. — Upper Panel: the SEDs for the different values of the high energy electron distribution 
index. The black dashed lines represent the synchrotron emission without self absorption. The 
black dotted line represents the synchrotron self adsorbed emission. The black solid lines corre- 
spond to SSC emission. The red and blue dotted lines correspond to the dusty torus and BBB 
emission respectively. The red and blue solid lines correspond to the ERC/DT and ERC/BLR 
emission respectively. Lower Panel: optical/7— ray flux correlations for ERC/BLR (blue solid cir- 
cles), ERC/DT (red solid circles), and SSC (black solid circles); dashed lines are the power law 
relations. 
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